NACA TN No. 1851 


. y3,t\2i/s:i/irsi 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 

No. 1851 


1 

► 


CRITICAL SHEAR STRESS OF INFINITELY LONG, SIMPLY 

SUPPORTED PLATE WITH TRANSVERSE STIFFENERS 

By Manuel Stein and Robert W. Fralich 

Langley Aeronautical Laboratory 
Langley Air Force Base, Va. 


I 



Washington 
April 1949 




APR II 1&4* 


8USINESS, SCIENCE 

& I echn^logy dept. 


NATIONAL ADVISORY COMMITTEE EOR AERONAUTICS 


TECHNICAL NOTE NO. I85I 


CRITICAL SHEAR STRESS OF INFINITELY LONG, SIMPLY 
SUPPORTED PLATE WITH TRANSVERSE STIFFENERS 
By Manuel Stein and Robert W. Fralich 

SUMMARY 


A theoretical solution is given for the critical shear stress of 
an infinitely long, simply supported, flat plate with identical, equally 
spaced, transverse stiffeners of zero torsional stiffness. Results are 
obtained by means of the Lagrangian multiplier method and are presented 
in the form of design charts. Experimental results are included and 
are found to be in good agreement with the theoretical results. 


INTRODUCTION 


The design of shear web beams and nonwrinkling skin surfaces 
requires a knowledge of the critical shear stress of stiffened plates. 
The purpose of the present paper is to give the theoretical critical 
shear stress of an infinitely long, simply supported, flat plate rein- 
forced with identical, equally spaced, transverse stiffeners. 

The results are found by means of the Lagrangian multiplier method. 
The stiffeners are assumed to have bending stiffness but no torsional 
stiffness and are assumed to be concentrated along transverse lines in 
the middle plane of the plate. The assumption that the stiffeners have 
no torsional stiffness applies with little error in the case of many 
open section stiffeners. The assumption that the stiffeners are con- 
centrated along transverse lines in the middle plane of the plate is 
applicable whenever the width of the attached flange is small in com- 
parison with the stiffener spacing. 

The theoretical analysis of the problem is given in the appendixes. 
For completeness, an energy solution for the plate with relatively weak 
stiffeners is given in appendix A. The solution for a plate with 
stiffeners of intermediate or higher bending stiffness is given in 
appendix B. The results are presented in the form of nondimensional 
curves which cover the complete range of stiffener stiffness and various 
stiffener spacings and in a table giving values from which the curves 
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■were drawn (table i). Experimental results are presented for 20 panels. 
Comparison of these results with the present theory indicates good 
agreement between theory and experiment. 


SYMBOLS 


T 


k 


s 


t 

b 

d 

b/d 


D 


E 

I 


P 

El 

Dd 


w 



critical shear stress 


critical sheaj>-stres3 coefficient 



thickness of the plate 
width of plate 
stiffener spacing 
panel aspect ratio 

flexural stiffness of the plate 


/ E_t 3 


\l2(l - p 2 ) 


Young’s modulus for plate 


Young’s modulus for stiffener 
effective moment of inertia of stiffener 
Poisson’s ratio for material 


ratio of stiffener stiffness to plate stiffness 


half wave length of buckles 
deflection of the plate 
deflection of the i’th stiffener 


x,y reference axes 

m, n, q^l 
r> 1 


integers 
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coefficients of deflection function 


undetermined Lagrangian multipliers 


V 


internal energy of bending of the plate 


V, 


s 


internal energy of tending of stiffeners 


T 


external work of the stresses 


BACKGROUND 


The problem of the buckling of stiffened plates in shear has been 
treated by many authors by the use of both theoretical and semi- 
empirical methods. In 1930 Schmieden (reference 1) solved the differ- 
ential equation for an infinitely long plate stiffened by closely 
spaced transverse stiffeners (equivalent to orthotropic plate) and 
found exact stability criterions for shear buckling of plates with 
simply supported edges and with clamped edges. By making certain 
simplifying modifications of the stability criterions, Schmieden 
obtained approximate values of the critical shear stresses. Later in 
1930 Seydel (reference 2) obtained exact solutions for infinitely long 
orthotropic plates with simply supported or clamped edges. With the 
use of the proper parameters Seydel* s results can be readily applied 
to plate— stiffener combinations. The values of the stresses obtained 
from Schmieden *s theory lie slightly below the exact values of Seydel. 
In 1947 T. K. Wang (reference 3 ) used the energy method to obtain an 
approximate solution for plate— stiffener combinations with simply 
supported edges. Wang’s results lie above the exact values of Seydel. 
All the foregoing solutions are applicable only to the case of weak 
stiffeners, where the stiffening effect of the stiffeners can be 
considered to be uniformly distributed over the plate. 

Solutions are also available for plates reinforced by rigid 
stiffeners. In 1936 Timoshenko (reference 4) treated the case of 
simply supported rectangular plates reinforced with one or two 
stiffeners. By means of the energy method Timoshenko found the 
stiffener flexural rigidity necessary to prevent buckling across 
stiffeners with the conservative assumption that the stiffeners act as 
simple supports. In 1948, Budiansky, Conner, and Stein (reference 5 ) 
found the critical shear stress for an infinitely long, clamped plate 
divided into square panels by nondeflecting intermediate supports which 
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correspond to rigid stiffeners. They also considered the case of a 
plate of infinite length and width having nondeflecting intermediate 
supports that form an array of square panels. 

Kuhn has written a number of papers on related subjects in which 
he presents semiempirical results for the critical shear stress of 
stiffened plates. (See, for example, reference 6.) 

The available theoretical solutions treat the relatively unim- 
portant case of weak- or closely spaced stiffeners and the case of 
rigid stiffeners that divide a plate into square panels. None of the 
theoretical solutions presents results for the practical range of 
intermediate stiffener stiffness and very little theory is presented 
for the practical range of spacing of rigid stiffeners. Also, it is 
felt that the semiempirical results for transverse stiffened plates 
cannot be extended to all stiffener spacings and stiffnesses without 
a sound theoretical basis. The theoretical results of the present 
paper cover the complete range of stiffener stiffness and the practical 
range of stiffener spacing. 


RESULTS AND DISCUSSION 


The critical shear stress for a plate — stiffener combination 
is given by the formula 


Curves are presented in figure 1 giving corresponding values of k Q 

TpT 

and the stiffness parameter — for simply supported, transversely 

Dd 

stiffened plates with panel aspect ratios of 1, 2, and 5- These 
results are replotted in logarithmic form in figure 2 for comparison 
with experimental results. 

The points of discontinuity of the slopes in the curves of figure 1 
represent changes in buckle patterns. The present results for an ortho- 
tropic plate agree with the exact results of reference 2. The deri- 
vation of the buckling criterion for an orthotropic plate (a plate 
stiffened by stiffeners of low bending stiffness) is given in appendix A. 
The derivation of the buckling criterion for plates stiffened by 
stiffeners of higher bending stiffness is given in appendix B. 
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In previous solutions, values of k g were found by using the 
orthotropic— plate curve and a cut-off at the value of k g for simply 

supported panels. (See fig. 1.) These figures show that the present 
solution yields values of k g that are considerably below those given 

by the orthotropic— plate curve in the intermediate range of stiffener 
stiffness. Also, the present solution for more rigid stiffeners yields 
a curve that is higher than the cut-off, which is obtained by assuming 
the stiffeners to have the effect of simple supports. Since the conti- 
nuity of the plate across the stiffeners of higher bending stiffness 
certainly adds a constraint to the plate, a higher buckling stress than 
that corresponding to a simply supported edge is obtained. 

In figure 2, experimental results are compared with the theoretical 
curves. These results are from two sources. The first set of experi- 
mental data is taken from MCA tests on shear webs of 24S-T aluminum 
alloy attached to torsion boxes. Drawings of a shear web and torsion 
box and the method of loading are given in reference 7- Buckling loads 
were obtained from the stiffener load-deflection curves which were 
taken from the original data. Each of the buckling loads given in the 
present paper is the average load at which the stiffeners start to 
deflect. The properties of the specimens and the buckling data are 
given in table II. 

The second set of experimental data is taken from MCA tests on 
thick web beams described in reference 9. The beams were made of 
24s— T aluminum alloy with heavy flanges and with joggled stiffeners 
riveted to the flanges. The open spaces in the joggles were filled 
with soft metal. A picture of a failed beam is shown in figure 3 . 

The load was applied at the center and the reactions were at the 
ends of the beams. Lateral deflections were prevented by lateral 
supports. The load, when strain was first observed in the stiffeners, 
was taken as the buckling load. The properties of the specimens and 
the buckling data are given in table III. 

The stiffener spacings for the test results are not the same as 
those for the theoretical results. All the test results fall in the 
expected regions among the theoretical curves. Only the group of test 

results for which = 2.4 fall in the range which serves to verify 

the present theory over previous theory which considered the orthotropic— 
plate curve to hold up to the cut-off at which the stiffeners are 
assumed to act as simple supports. The other groups of test results 
agree with the present theory, but they do not cover the range in which 
an appreciable difference exists between the present theory and previous 
theory. More experimental results are required to confirm the present 
theory fully. 
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CONCLUDING REMARKS 


Charts are presented from which the theoretical critical shear 
stresses can he obtained for infinitely long, simply supported plates 
stiffened with identical, equally spaced, transverse stiffeners of zero 
torsional stiffness. The theoretical results are based on the Lagrangian 
multiplier method. Previous theory considered the orthotropic curve 
to hold up to a cut-off value corresponding to the stiffener stiffness 
at which the buckling load was equal to the buckling load of a simply 
supported panel the size of each bay. Comparison of the present theory 
and previous theory shows that previous theory gives unconservative 
results for stiffeners of intermediate stiffness and conservative 
results for stiffeners of high stiffness. Test results of 20 panels 
are presented which are in good agreement with the present theory. For 
a conclusive check additional test results are required. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., January 28, 19^9 


NACA TN No. I 85 I 


7 


APEENDIX A 


THEORETICAL SOLUTION OF CRITICAL SHEAR STRESS OF ELATES* 
WITH TRANSVERSE STIFFENERS OF LOW BENDING STIFFNESS 


If the stiffener "bending stiffness is low and the stiffeners are 
fairly closely spaced, the buckle pattern may be considered independent 
of the stiffener spacing, and the plate stiffener combination can then 
be analyzed as a plate with different bending properties in each 
direction, that is, an orthotropic plate. In this appendix buckling in 
shear of an orthotropic plate is analyzed by means of the energy method. 

The buckling configuration of the plate shown in figure 4 is 
represented by the trigonometric series 


which satisfies the boundary conditions of simple support term by term. 
The internal bending energy of the plate V, the internal bending energy 
of the stiffeners V , and the external work of the shear stresses T 

are given by the expressions 


00 


00 



y Jy • • • 



(Al) 




(A2) 



J 
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Substitution of the expansion for w (equation (Al)) into these 
energy integrals gives 



T = 2rtrt 



n— 1»3 3 • • • 



a„b„ 

q n 



Then 


(V + V s - T) 


8b3_ 


00 




n^ — 

Dd 


00 

n— 1^3^ • • • 




n=l,3> ••• q=2j^j 


a q^n 


nq 


(A3) 


n c - q 
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where 


k = 
s 


t tb £ 
D * 2 


According to the energy method the potential energy (y + y - T) 
mnat he minimized with respect to the unknown coefficients a n S and b 
By minimizing (v + v g - T) with respect to the coefficients a and b. 
the following set of equations is obtained: n 



81*3 

nX 


OO 

n 

<j— • • • 



. . . 


nq 


= 0 


11 (l 2 - n2) 

(n=2,k,6 , . . . ) 


nq 

a = o 

3 (- 2 -* 2 ) 

3j 5* • • • ) 


(a4) 


(A5) 


The coefficients a n can be found in terms of b r from 
equation (Ah). Substitution of the resulting expression for a in 
equation (A5) results in the following equations: 
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A solution to equations (A6) exists if the following stability determi- 
nant vanishes: 


C n 

°13 

c 15 

c 3i 

c 33 

c 35 

c , 

c 

C 

51 

53 

55 

• 

• 

• 


(AT) 


where 



(where n ^ r) 


A solution including all the a * s and b_ can be obtained by 

n 1 

setting equal to zero the first approximation of the determinant 
equation (A7) • 


C 


11 


= 0 
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Similarly the second approximation includes a] ~l the a^’ s, 
and b^ 


b 


1 ' 


C 


11 



= 0 


Higher approximations are found in a similar manner. A second approxi- 
mation was found to give satisfactory results. For a given approxi- 
mation it is necessary to try values of b/X and find the corre- 
sponding values of k g until a minimum value of k g with respect 

to b/X is found for each 51. The results are given in table I and 


in figure 1 
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APPENDIX B 

THEORETICAL SOLUTION OF CRITICAL SHEAR STRESS OF PLATES 
WITH TRANSVERSE STIFFENERS CF HIGHER SENDING STIFFNESS 


In appendix A a theoretical solution for a plate stiffened by 
stiffeners of low bending stiffness is presented where the buckle 
pattern was taken as sinusoidal in the longitudinal direction. The 
buckle pattern of plates with stiffeners of higher bending stiffness 
is. no longer sinusoidal in the longitudinal direction. It is then 
necessary to consider deflection functions which are either symmetric 
or antisymmetric about the midpoint of each bay and are periodic over 
an integral number of bays. The critical shear stress of plates with 
transverse stiffeners of higher bending stiffness is analyzed by means 
of the Lagrangian multiplier method. 

Deflection functions. — The correct buckle configuration for any 
given plate-stiffener combination is that which corresponds to the 
lowest buckling load. Several types of configurations are investi- 
gated. These buckling configurations are represented by the following 
two-dimensional trigonometric series (the coordinates are given in 
fig. 4). Symmetric buckling, periodic over each bay: 


w 


-XI X 


m=2,4, ... n=2,k, . 


. mrtx . nrty 
a sin -3— sin — — 
mn d -h 


x: x 


. mrtx . nrty 

b^ cos — — sin 
mn d b 


m— 0, 2, . . . n-1, 3 y ' 


Antisymmetric buckling, periodic over each bay: 


v = > 2 


. mrtx , nrty 
a sin — — sin —r— 
mn d b 


m=2,4, ... n=l,3,, 


(Bla) 


XZ X 


, mnx . n«y 

b cos sin — — 

mn fl -h 


m=0,2 , ... n=2,4,, 


(Bib) 
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Symmetric buckling, periodic over two bays: 


as 00 



m— 1 , 3 i • • • H— 1 , 3 , • • • 


a 

ran 


. mrtx . nity 

am 3m — — 

d b 


2Z 2 


m=l,3,... n=2,4. 


, mrtx . nity 

b mn cos — sin 


(Blc) 


Antisymmetric buckling* periodic over two "bays: 


w = 


-51 2 


o . mrtx . nity 
a mn sin — sin 


m=l,3.>... n=2,4, , 


5: 2 


m— 1 y 3 f • • • n— 1 j3i • • • 


mrtx nrty 
b cos sin 

mn d b 


(Bid) 


Symmetric buckling, one .bay; antisymmetric buckling, next bay; periodic 
over four bays: 


w = )> "> 


a mn 


m=l,3, ... n=l,3 , ... 


m-1 

sin 22£ + (_!) 2 coa mtx 
2d ' ; 2d 


sin 22Z 
b 


51 5 


mn 


m=l,3,... n=2,4, , 


m-1 

mrtx / -1 \ 2 m» tx 

sin — (—1) cos 

2d 2d 


sin 


n«y 


(Ble) 
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Careful study has shewn that other "buckle patterns would require higher 
buckling loads and that only the five buckle patterns given need be 
considered. 

These deflection functions all satisfy term by term the conditions 
of simply supported edges at y = 0, b and continuity of the plate 
across the stiffeners at 1 = 0, d, 2d,.... The condition that 
stiffener deflection equal plate deflection at the stiffeners is 
introduced by means of Lagrangian multipliers. 

The deflection functions (Bid) and (Ble) are found to be the 
governing ones for the aspect ratios investigated; the others lead 
to unconservative solutions. Buckling criterions for the critical 
shear stress are derived for the deflection functions (Bid) and (Ble). 

Antisymmetric buckling, periodic over two bays.— The deflection 
of the plate is given by equation (Bid) as 


00 CO 



nt=l,3, • • • n— 2,4, . . . 


00 00 
m— 1, 3 , . o . n — 1 , 3 ) • • • 


mrti nrty 
b cos sin 

mn d b 


The deflection of the i*th stiffener is taken as 


00 

( w b)i ‘ 2 . ^>1 sln T 

n=l,3, • • • 


(B2) 
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where, since the interval to he considered includes two stiffeners, i = 1 
and 2. The boundary conditions that stiffener deflection equal plate 
deflection are 


v(id,y) - (v^ = 0 


(i=l,2) 


or upon substitution. 


> b mn + ^l = 0 (n=l,3,...) 

3°° b mn -A n2 = 0 (n=l,3,...) 

m = lj 3 ) • • • 


These equations show that A^ = -A^. If A^ is redefined as A n 
the boundary conditions become 


1> Vi + \ = 0 ( n =lj3j . . . ) (B3) 

m=l,3,... 

These boundary conditions will be satisfied in the energy expression by 
means of Lagrangian multipliers. 

The internal bending energy of the plate V, the internal bending 
energy of the stiffeners V s , and the external work of the shear 

stresses T are given by the expressions 
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Substitution of the deflection functions of the plate and stiffeners 
into these energy integrals gives 


V = 


Ddit 

4b 3 


x 2 


HI— 1 y 3 9 • • • n-2 y ^ y • « . 


n 2/ m 2b^ p' 2 

Ann l m ^7 




^ 2 / o -u 2 0 \2 

b / t . -2 


mn 


m — + n 


hi— lj 3 s • • • n-1, 3, . . . 


V a = 

8 2b 3 


1 CD 

El* ^ a 2 4 

Ai n 


u— l^ 3,... 


)> (B5) 


T = ^Ttrt 


X X 2 


a b __JE3, 

mn mq 


HI— 1 y 3 y • • • HI— 2 ; ^t- j • q— 1 ^ 3 y © • < 


(q 2 - n 2 .) 


J 
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The energy method requires that the potential energy (V + V g - T) 

he minimized with respect to the a’s, b*s, and A*s. Since the a’s, 
h’s, and A’s are, however, hound hy equations (B3), the minimization 
is performed hy the Lagrangian multiplier method hy minimizing the 
following function F with respect to the a’s, h*s, and A’s: 


F = 


V + V 3 - T 
/ tc^D d \~ 

\ 4h3 ) 


00 / 00 

^n ( ^ ^mn + Aa 

Vm=l,3,..- 


(B6) 


n— 1,3^ ... 


where the 7 ’s are the Lagrangian multipliers. When this minimization 
is performed, the following set of equations is obtained: 


dF 

Ka 


2 b £ 


= 0 = 2a mn m T + n 


\ 2 16 k, 


—1 £ 
jc d y 


q-i^3^ • • < 


8 f 


8 b 


mn 


, 0 -k^ 2 \ 2 l6k, 

0 = 2 b m 2 S- + n + 
mn \ ^2 


s h 

* d 


I 


q— 2,4, • , 


h 7 — SIS — ^ 

“4 („a _>) 


(m— 3 9 • • • ) 

(n«2, k,. 9 .) 


mnq 

°“ q * r ° 

(m=l, 3 , . . . ) 

(n=l, 3 , • • • ) 


(B7) 


8 f _ 4ei 4 
- — = 0 = n A_ + 7 

8 ^ Dd " n 


( n= l,3, • . . ) 


J 
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When the equations (B7) are combined., the following equations are 
obtained : 



Equations (B 8 a) written in matrix form are 







C ml C ml3 C ml5 


b ml 


2 A, 

C m31 c m3 C m35 * * * 



= El 

I 62 A 3 

C m51 C m53 C m5 

« • • 


b m5 

• 

M 

1250 ^ 

• 

• • • 

• • • 


• 

• 


• 

• 


(m— 1 , 3 , . . . ) 


where 



(B 8 a) 


(B 8 b) 
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A solution including all the und b^ * a can 1© obtained 

by the first approximation of the matrix equation (B8b) 

C ml^ml ~ ^ ^1 (m=l,3j • • • ) (B9) 

» 


Substitution of b^ from equation (B9) into the boundary equation (B3) 
yields 



(BIO) 


The following stability criterion is obtained by setting equal to 
zero the coefficient of Aj_: 



m= l*3j • • • 



(Bll) 


Similarly, the second approximation includes all the a »s, b »s 

mn 3 ml J 

and b f s. Tvo simultaneous equations result from which b and b 

ml m3 

can be found. Substitution of these values into the boundary equation (B3) 
yields two linear homogeneous equations in and A^. If the determi- 

nant of the coefficients of these two equations is set equal to zero, the 
following stability criterion is obtained: 
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"m!3 


C ml C m3 “ C ml3 4 


(B12) 


Higher approximations are found, in a similar manner. A second 
approximation vas found to give satisfactory results. For each of 
these approximations, it is necessary to find the lowest value 

El 

of lr for each value of — . The results are given in table I and 
3 Dd 

in figure 1. 

Buckling periodic over four hays.- The deflection of the plate is 
given by equation (Ble) as 


c» 


03 


-XL 5 


mn 


ni—1 j 3 • • ' 


m-1 

nan , . 2 hbcc 

sin + (—1) cos 

2d 2d 


sin 


iucx 


00 




mn 


nt-l ,3 , o • . n-2 ,4, • • • 


mrtx , . \ 2 mux 

sin (-1) cos 

2d 2d 


sin 


nxy 


The deflection of the i’th stiffener is taken as 


n*l, 2, . . • 


sin 22Z. 
b 


(B13) 


MCA TN No. 1851 


21 


where i = 
stiffeners. 


hecome 


2, 3 ) and. 4, since the interval considered includes four 
The boundary conditions 


v(id,y) - (v s )i = 0 


(i-1,2,3,4) 


°q m-1 

a mn “ ^nl = 0 ( n=1 ,3, ...) 

m. = lj 3 y ... 

°° m-1 

^ ^mn = ® (n=2,4 x . . . ) 

m?=l y • • • 


> a nm (- 1 ) 2 + \,2 - 0 


m-1 

2~ 


nfc=lj3 , . . . 


(n— 1^ 3* • • • ) 


^ ^mn (“•!) ^ *” Aq2 ” ® 


m-1 

2 “ 


m^l ^ • • • 


(n-2,4, • • • ) 


CD 

nk-lj 3* < 


m-1 

^ + ^3 = ° 


(n=l,3,...) 


]>Z v <- i) 2 + ^3 = 0 


m-1 

2 


m-1 y 3 y • • • 


(n-2j4, • • . ) 
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00 m — 1 

^ “am ( -1 ) “ = 0 

m=l y 3 y • • • 

00 gfc— 1 

>— v, (- 1 * 2 + ^ ° 

EQr=l y 3 y • • • 

These equations show that 


A nl ^h2 ^3 " 

(n-1,3 

= Ai2 = " ^3 = ~ 

(n=2,4 


(n=l,3,...) 


(n— 2,4> • . . ) 


If la redefined as the boundary conditions become 


5 


mr-1 

a nm (- 1 ) 2 ~\ = 0 


(n— 1, 3* • • •) 


> (B14) 


00 nt — 1 

y~ v (- 1 ) 2 - ^ - ° 


(n-2,4, . . • ) 




These boundary conditions vill be satisfied in the energy expression by 
neans of Lagrangian multipliers. 

The energy integrals are the same as the energy integrals (B4 ) , 
except that in the present problem the upper limit of integration 2d 
is replaced by 4d and the upper limit of the summation 2 is 
replaced by 4. 
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The deflection functions of the plate (equation (Ble)) and stiffeners 
(equation (BI 3 )) are substituted into these energy integrals and result 
in the following expressions: 


V = 


Debt ' 

b3 


00 ~ 


m- 1 , 3, • • • ti—1,3j • • • 


2/m 2 b 2 2> 

a mn U d 2 + n J 


x: x 


21 — 1^3 j • • • n — ... 


2A l2 b 2 2\ 2 

tF + n ) 


V B - 


EIjT 

b3 




n=l, 2 , . . . 


T = drtn 


H X 2 


m— 1,3* • • • ft~l,3j • • • m • • 


m -1 

a b (- 1 ) 2 —S3 

■"■3 r n S _ 2\ 


The minimization of (V + V 8 — T) is performed by the Lagrangian 

multiplier method by minimizing the following function F with respect 
to the a*s, b’s, and A*s. 


F = 


V + V a - T 
rt^Dd 

b3 


J 


'n 


n=l, 3 ,... 


2 


m -1 

%n (-1) 2 -Aa 


m- 1 , 3> • • 


2 : 


n 


n=2,4, . . . 


00 

21—1 y 3 y • • • 


m -1 


V (- 1 ) - ^ 


(B15) 


where the 7*3 are the Lagrangian multipliers* 
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When the minimization is performed and the resulting equations 
are combined, the following set of equations is obtained: 


■\m a - 


mn m 




(j— 2^4^ < 


Itt-1 

(m=lj 3* • • • ) 

(n=l^ 3* • • • ) 


■^mn^mn ^m ^ 


nq 


Q-l-s 3 * 1 




H)r— 1 

+ 2 — nV (-1) 2 = 0 


M 


(m— 1, 3, . . . ) 
(n=2,4, ...) 


> 


y 


where 


A 


mn 



r m 


8k "b 
s 

rt& 


mr-1 

m(-l)"^~ 


% 


(Bl6a) 
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Equations (Bl6a) in matrix form are 


Anl 


- r 

3 m 


0 

iLr 

15 m 

0 


- r m 0 

3 1 

V -f r m 

- - r a 

5 m m3 


15 


m 


Ar 

35 


m 


-A r 

35 m 


0 

_ 10 r 

" 21 m 
0 


12 r 


m 


- r 


12 r 

•j m 


An.^ 


10 r 
” 2.1 m 

0 


20 p 

7 m 


- r 

3 m 


— p 


m 


V 5 u F m 


m 


11 m 


■\i6 


*mL 


D m2 


a. -v 

m3 


3 m4 


a 


m5 


m6 


m+1 

(-U~ H 

Dd 


2 Al 


32 Ag 


162 


*3 


512 ^ 

1250 A^ 

2592 A6 


(Bl6b) 


(m— 1, 3, . . . ) 


A first approximation of k a is found by considering all the a^’s 


and b^’s in equation (Bl6b). 


X 


m+1 


VlVl + f r ,V - (-D 2 2 § A 1 (®»U3. • • ■ ) 


m+1 

O 'p t?t 

J ^m a ral + AnP^mP = ( — 1) 32 — Ag (m=l, 3* • • • ) 




(B17) 


Substitution of a^ and b^ from equations (B17) into the boundary 
equations (Blk) yields 
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1+2 .El 
1 m=l, 3, „ . . 9 ^ 


1 «■ 


1*1 


!> 


hi 


m=l_, 3, . . . Vi/W 


4 „ 2 


i*2 = 0 


- 2. r. 


m 


(B18) 




m 


m— 1^ 3.> • • 


. AmlAm2 9 *i& 


1*1 


ml 


un— 1,3, • • 


*1 2 

g m 


F2 


= 0 


V\2 - - ^ 32 


If the determinant of the coefficients of the linear homogeneous 
equations (Bl8) is set equal to zero, the following stability criterion 
is obtained; 



Similarly, from th9 second approximation, including all the a ml' 

b and a - terms, the following stability criterion is obtained: 
Hki -7 m3 


(B19) 


36 2 

a a — — r 

m2 m3 29 m 1 


B, 


2 r mV 3 


m=l,3, . . . 


m 


_ El 3 0 B 

2 *— m-l, 3, . . . in 


where 


^nAn3 


Dd 



Bm 

X B 

m=l , 3 , . . . 

in 

m^l, 3 , • . . 

), “ 

r 2 

m 

6 x 00 

5 ]> 

5> 

m-1, 3* • • • 

IDr-1, • • 


+ _i. 


32 S 
Dd 


^m A ml 

B m 


4 ^ 

“ ^ B* 


m— 1 j 3 y • • * 


is 


r a 

m ml 

^m 




3^ • • • 

k 2 

^ml^in2 9 ^in. 
B 

m 


162 Si 

Dd 


= 0 (B 20 ) 


B = A A A , - — r 2 A -- r \ 
m ml m 2 m 3 25 m 9 m m3 


Xn a similar manner., a stability criterion is obtained from a third approximation, which includes 
all the a rn1 , b^, a^, and b^ terms. A second approximation was found to give satisfactory results 

for most cases. For certain cases noted in table I far a panel aspect ratio of 5, however, third 
approximations should, be used. The terms used in these approximations are given in table I. For each 

of these approximations, it is necessary to find the lowest value of k for each value of — . The 

3 Dd 

results are given in table I and in figure 1. 
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TABLE I. - SHEAR— STRESS COEFFICIENTS FOR STIFFENED PLATES 
WITH PANEL ASPECT RATIOS OF ONE, TWO, AND FIVE 


Plates with stiffeners of low 
bending stiffness 8. 

Plates with stiffeners of higher 
bending stiffness 

El 

Dd 

k 

s 

Aspect 

ratio 

Antisymmetric 
buckling 
periodic over 
two bays 

Buckling 

periodic 

over 

four bays 

11 

Dd 

k 

s 

El 

Dd 

k 

s 

0 

2 

5 

20 

50 

100 

200 

5-3^ 

10.34 

16.07 

37-14 

68.99 

112.2 

184.6 

1 

0 

2.91 

7.78 

22.29 

00 

5-53 

7.85 

9.82 

H.78 

13.86 

0 

7.09 

19.03 

00 

6.08 

10.0 

10.5 

10.86 

2 

0 

3-35 

14.50 

22.99 

33.11 

45.77 

61.97 

82.92 

112.3 

605 

00 

9.65 

12.0 

16.0 

18.0 

20.0 

22.0 

24.0 

26.0 
28.0 
35.0 
37.05 

0 

5.475 

11.93 

26.37 

36.29 

68.92 

145.4 

625 

00 

5.54 

15.0 

20.0 

23.0 

24.5 

26.0 

27.0 

28.0 

28.2 

5 

0 

18.02 

90.99 

176.8 

444.7 

704.4 

42.5 

70 

90 

100 

120 

140 

0 

49.19 

112.8 

220 

00 

13.37 

60 

b ioo 

c i4o 

c i43 


independent of aspect ratio. 

All the ami, bj^, a^, and b^ coefficients used. 
C A11 the ajj^, hjjjj^, a^, and b^ coefficients used. 
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TABLE II.- EXPERIMENTAL BUCKLING DATA OF SHEAE WEBS 
WITH UPRIGHTS NOT CONNECTED TO THE FLANGES 


Specimen 

(a) 

d 

(in.) 

b 

(in. ) 

t 

(in.) 

Uprights 
(nominal size) 
(in.) 

El 

Dd 

T 

(ksi) 

k s 

2-D-O 

5-0 

23-5 

0.0397 

1/2 x l/2 x l/l6 

221 

2.66 

101 

3-D-o 

5-0 

23.5 

.0394 

3/4 X 3/4 X 1/16 

680 

3.08 

116.5 

4-D-O 

5-0 

23.5 

.0405 

3/4 X 3/4 X 3/32 

946 

3.295 

U7.5 

5-D-0 

10.0 

23.5 

.0404 

1/2 X 1/2 X 1/16 

98.3 

1.21 

43.3 

6-d-o 

10.0 

23.5 

.o4o8 

3/4 X 3/4 X 1/16 

306 

1.54 

54.2 

T-d-o 

10.0 

23.5 

.0410 

3/4 x 3/4 x 3/32 

456 

1.47 

51.3 

8-S-O 

5-0 

23.5 

.0394 

1/2 x 1/2 x 0.064 

95-8 

2.895 

109 

9-s-o 

5-0 

23.5 

.0399 

3/4 x 3/4 x 3/32 

456 

3.01 

111 

10-S-0 

10.0 

23.5 

.0410 

1/2 X 1/2 X 1/16 

41.4 

.82 

28.6 

11-S-O 

10.0 

23.5 

.0398 

3/4 X 3/4 X 1/16 

151.5 

1.357 

50.1 

12-S-O 

10.0 

23-5 

.0405 

3/4 X 3/4 X 3/32 

217 

1.41 

50.3 


a S, stiffeners on one side of plate. 

D, stiffeners on Both sides of plate. 
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TABLE III . - EXPERIMENTAL BUCKLING DATA OF THICK WEB BEAMS 
WITH UPRIGHTS CONNECTED TO THE FLANGES 


Specimen 

(a) 

b 

(in.) 

d 

(in.) 

t 

(in.) 

Uprights 
(nominal size) 
(in.) 

El 

Dd 

T 

(Esi) 

k 

s 

V- 12 - 7 S 

9.88 

7-00 

0.1005 

li X li X i 
8 8 8 

91.0 

15.5 

15.4 

V- 12 - 8 S 

9.88 

7.00 

.1044 

3/4 x 3/4 x 1/8 

25-8 

15.4 

14.15 

V- 12 - 9 D 

9.13 

7.00 

.1025 

5/8 x 5/8 x 1/8 

40.4 

16.8 

13.65 

V— 12 — 10 S 

9.88 

7.00 

.1043 

5/8 x 5/8 x 1/8 

14.5 

I6.3 

15,0 

V— 12 — 11 D 

9.13 

7.00 

.1025 

5/8 x 5/8 x 3/32 

30.3 

17.2 

14.0 

V— 12 — 12 S 

9.88 

7.00 

.0987 

1/2 x 1/2 x l/l6 

4.1 

12.3 

12.7 

V- 12 - 13 D 

9.13 

7.00 

.1000 

1/2 X 1/2 x 1/16 

11.3 

13.1 

11.15 

V— 12 — 14 S 

9.88 

7.00 

.1007 

5/8 x 5/8 x 3/32 

11.2 

13.2 

13.1 

V- 12 - 15 D 
] 

9.13 

7.00 

.1057 

5/8 x 5/8 x 1/16 

18.8 

15.7 

12.0 


S, stiffeners on one side of plate. 

D, stiffeners on Both sides of plate. 
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(t>) Panel aspect ratio — = 2. 
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Figure 1.- Continued. u> 
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(c) Panel aspect ratio 



Figure 1.- Concluded. 


NACA TN No. 1851 


10 


Figure 2 
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Figure 3 .- Thick web beam (with uprights connected to flanges) after failure. 
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Figure 4 .- Infinitely long, simply supported plate, with transverse 

stiffeners, under shear. 


